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Transverse mass and rapidity distributions for charged pions, charged kaons, protons and antipro- 
tons are reported for ^/snn = 200 GeV pp and Au+Au collisions at RHIC. The transverse mass 
distributions are rapidity independent within \y\ < 0.5, consistent with a boost-invariant system in 
this rapidity interval. Spectral shapes and relative particle yields are similar in pp and peripheral 
Au+Au collisions and change smoothly to central Au+Au collisions. No centrality dependence was 
observed in the kaon and antiproton production rates relative to the pion production rate from 
medium-central to central collisions. Chemical and kinetic equilibrium model fits to our data reveal 
strong radial flow and relatively long duration from chemical to kinetic freeze-out in central Au+Au 
collisions. The chemical freeze-out temperature appears to be independent of initial conditions at 
RHIC energies. 

PACS numbers: PACS number(s):25.75.Dw 

Quantum chromodynamics predicts the existence of achieved in relativistic heavy ion collisions . Signals of 
a new form of matter, the quark-gluon plasma (QGP), QGP may remain in the bulk properties of the collision, 
at extreme conditions of high energy density, possibly and simultaneous observations of multiple QGP signals 
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in the final state would serve as a strong evidence of QGP 
formation. These bulk properties include strangeness and 
baryon production rates and collective transverse radial 
flow. These can be studied via particle spectra. 

In this letter we report results on charged pion (w ), 
charged kaon (K ), proton (p) and antiproton (p) pro- 
duction from pp and Au+Au collisions at RHIC by the 
STAR experiment at the nucleon-nucleon center-of-mass 
energy of ^/snn = 200 GeV. In some models it is argued 
that particle multiplicity density per transverse area of 
interaction measures the initial gluon density 0], parti- 
cle ratios measure the chemical freeze-out conditions 0, 
and transverse momentum spectra measure the kinetic 
freeze-out conditions Q. We study these properties at 
mid-rapidity as a function of centrality. The rapidity de- 
pendences of particle production and spectra shape are 
also investigated. 

Charged particles are detected in the STAR Time Pro- 
jection Chamber (TPC) 5]. The TPC is surrounded by a 
solenoidal magnet providing a uniform magnetic field of 
0.5 T along the beam line. Zero degree calorimeters and 
beam-beam counters provide a minimum bias trigger 
for Au+Au and pp collisions, respectively. Events with a 
primary vertex within ±25 cm of the geometric center of 
the TPC along the beam axis are accepted. For this anal- 
ysis, about 2.0 x 10 6 Au+Au and about 2.5 x 10 6 pp mini- 
mum bias accepted events are used. Only primary tracks 
- tracks pointing to the primary vertex within 3 cm - are 
selected. The Au+Au events are divided into 9 centrality 
classes based on measured charged particle multiplicity 
within pseudo-rapidity |7/| < 0.5. These classes consist, 
from central to peripheral, of 0-5%, 5-10%, 10-20%, 20- 
30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80% of 
the geometrical cross-section. 

Particle identification is accomplished by measuring 
the ionization energy loss dE/dx. The mean (dE/dx) 
is determined from 70% of the samples with the low- 
est dE/dx along a track. To insure good momentum 
and (dE/dx) resolution tracks are required to have at 
least 25 out of the maximum 45 hits in the TPC. The 
(dE/dx) resolution varies between 6% and 10% from pp 
to central Au+Au events. The reconstructed momenta 
are corrected for most likely energy loss in the detector. 
The correction is negligible for tt , under 2% for and 
under 5% for p and p in the covered momentum ranges. 
The momentum resolution was estimated to be about 2% 
at p± = 0.5 GeV/c. Uncorrected particle yields are ex- 
tracted from (dE/dx) distributions for ea,chp±, rapidity 
and centrality bin [tI IsL IqL fTo| . 

Corrections are applied to account for tracking inef- 
ficiency, detector acceptance, hadronic interactions, and 
particle decays. The total reconstruction efficiencies are 
obtained from embedding Monte Carlo (MC) tracks into 
real events at the raw data level and subsequently recon- 
structing these events. The propagation of single tracks 
is calculated using GEANT, a detailed description of the 



STAR geometry, and a realistic simulation of the TPC 
response 0, |l 03 • The efficiencies for tt* are 50-70% 
and 80-90% in the covered p x for the 0-5% and 70-80% 
events, respectively. The corresponding efficiencies for 
K± are 40-70% and 20-50% and for p and p 70-75% and 
75-80%. Background protons knocked out from the de- 
tector material are subtracted. This background is 50- 
60% at p± = 0.4 GeV/c and becomes less than 5% at 
1.0 GeV/c @|. 

Corrections for the pp data are similar to those for the 
70-80% Au+Au events. Additional corrections are ap- 
plied for primary vertex reconstruction inefficiency and 
fake events (events with mis-reconstructed vertex due 
to pile- up background). These corrections are obtained 
by embedding HIJING Tl| events into events that had 
been triggered on empty bunches, and reconstructing the 
combined events. The vertex reconstruction inefficiency 
strongly decreases with increasing event multiplicity re- 
sulting in approximately 14% of events being missed, 
over 80% of which have fewer than three tracks in the 
TPC. About 12% of pp events are fake events with re- 
constructed multiplicity about half of that of real events 
due to time distortion in the pile-up background, result- 
ing in an overall correction of 6-8% in the covered p± 
range. 

The pion spectra are further corrected for weak de- 
cay products, muon contamination and background pi- 
ons produced in the detector material. The weak decay 
correction is approximately 12% and was estimated from 
the measured K® and A distributions 0, extrapo- 
lated to our energy. Because weak decay (anti)protons 
carry most of the parent momentum, their tracks behave 
as those originating from the primary vertex, resulting 
in the same reconstruction efficiency for weak decay and 
primary (anti)protons over the measured p± range. 

The inclusive (anti)protons closely reflect total 
(anti)baryon production |j],|£|]. Therefore we present in- 
clusive proton and antiproton distributions that are not 
corrected for weak decays. Based on the measured A dis- 
tribution we estimate that about 40% of the mea- 
sured protons are from weak decays, and the measured 
inclusive (p±_) are similar to those of primary protons. 

The point to point systematic uncertainties on the 
spectra are estimated by varying event and track selec- 
tion and analysis cuts and by assessing sample purity 
from the dE/dx measurement. The estimated uncertain- 
ties are less than 4% for ir^, p and p. Those for 
are less than 12% for p± bins with significant overlap in 
dE/dx with or 7r ± , and less than 4% for other bins. 
An additional systematic error on the proton spectra due 
to background subtraction is estimated to be 5% at low 
p± and negligible at high p± ||. A correlated system- 
atic uncertainty of 5% is estimated for all spectra and is 
dominated by uncertainties in the MC determination of 
reconstruction efficiencies. 

Figure 1 shows transverse mass (m± = \/p\ + m 2 ) 
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m T - m „ ( GeV/c 2 ) m T - m K ( GeV/c 2 ) m T - m p ( GeV/c 2 ) 

FIG. 1: Invariant yield as function of transverse mass for ir , K , and inclusive p and p at mid-rapidity (\y\ < 0.1) for pp 
(bottom) and Au+Au events from 70-80% (second bottom) to the 0-5% centrality bin (top). Statistical and point-to-point 
systematic errors have been added in quadrature. Additional correlated systematic error due to uncertainty in the normalization 
is estimated to be 5%. Open circles are for positive particles (all proton spectra are scaled by 0.8), and closed triangles are for 
negative particles. The curves shown (Bose-Einstein fits for n~ and blast-wave model fits for K~ and p) are explained in the 
text. 



spectra for 7r ± , , p and p for pp and all centrality bins 
of Au+Au data within \y\ < 0.1. For clarity, proton spec- 
tra are scaled by 0.8. Particle and anti-particle spectra 
shapes are similar for each centrality bin. While the 
spectra shapes are similar for pp and Au+Au, K ± , p and 
p spectra show a progressive flattening from pp to cen- 
tral Au+Au events. Our pp results are consistent with 
previous measurements at similar multiplicities |l3l |. 

The blast-wave model — a hydrodynamically moti- 
vated model with a kinetic freeze-out temperature T^ in 
and a transverse flow velocity field /3 Q — can simul- 
taneously fit the K ± , p and p spectra and the h\gh-p± 
part (p± > 0.50 GeV/c) of the ir ± spectra. We used 
a velocity profile of (3 — [3 s (r/R) n , where r < R (the 
term r/R accounts for the change in the velocity as a 
function of radial distance), (3 S is the surface velocity, 
and n is treated as a free parameter. The value of n 
ranges from 1.50 ± 0.29 in peripheral to 0.82 ± 0.02 in 
central events. The fit results are superimposed in Fig. 
l(b,c). The obtained fit parameters for the 0-5% Au+Au 
events are T kin = 89+10 MeV and (/3) = 0.59 ± 0.05, 
P s = 0.84 ± 0.07, and are similar to the 130 GeV results 
reported in |9j, [T4| . The systematic uncertainties in the 
fit parameters are estimated by excluding the kaon or the 
(anti) proton spectra from the fit. 

Recent attempts to fit the measured RHIC spectra 
with a single (chemical and kinetic) freeze-out temper- 
ature claim this is possible if all the resonance and weak 



decay feed-downs are taken into account 15]. Our MC 
study of that scenario shows significantly higher x 2 /NDF 
compared to our blast-wave fits. 

The low-pj^ part of the pion spectrum deviates from 
the blast-wave model description, possibly due to large 
contributions from resonances at low p± . We fit the pion 
spectra to a Bose-Einstein distribution (oc l/(exp^ — 
1)), the results of which are superimposed in Fig. 1(a). 
The yields outside the measured p± region are extrapo- 
lated using the blast- wave model for K ± , p and p and 
the Bose-Einstein distribution for tt ± . The uncertainties 
on these extrapolations are estimated by comparing to 
results using other functional forms. The estimated ex- 
trapolation uncertainties in the (p_i_) and total yield are 
5% for tt^~ and 5 to 10% for K^, p and p (varying from 
pp to central Au+Au collisions). For the 0-5% Au+Au 
collisions, the integrated yields are dN/dy = 322 ± 32 for 
7T+, 327 ± 33 for n~, 51.3 ± 7.7 for K+, 49.5 ± 7.4 for 
K-, 34.7 ± 6.2 for p and 26.7 ± 4.0 for p. The obtained 
p/p ratio for the 0-5% Au+Au collisions is 0.77 ± 0.05, 
indicating a nearly net-baryon free mid-rapidity region 
at this RHIC energy. 

We extract the fiducial dN/dy by summing up the 
yields within the p± range of 0.20-0.70 GeV/c for 7r~, 
0.25-0.60 GeV/c for K~, and 0.50-1.05 GeV/c for p. 
Figure 2 depicts the rapidity dependence of the fiducial 
dN/dy and extrapolated (pj_) for the 0-5% and 70-80% 
Au+Au events. We do not observe changes in cither 
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shape or yield for any particle species within \y\ < 0.5. 
The pp data and all other centrality bins of the Au+Au 
data exhibit the same behavior. Such an absence of ra- 
pidity dependence of particle spectra was also observed 
for 7T^, v and p at ^/snn = 130 GeV Au+Au colli- 
sions [3, llj ■ This uniformity indicates the development of 
a boost-invariant region within the measured kinematic 
ranges. 
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FIG. 2: Rapidity distributions of the fiducial yields and inte- 
grated (px) for the 0-5% and 70-80% Au+Au collisions. Pion 
yields are scaled down by a factor of 5. Errors shown are 
those propagated from Fig.l. Systematic errors on the fidu- 
cial yields are 5%; those on (p±) are 5% for pions and 5-10% 
for kaons and antiprotons. 

The centrality dependence of the extracted (p±) within 
\y\ < 0.1 is shown in Fig. 3(a). A smooth changeover 
from pp to peripheral Au+Au collisions is observed for 
all particle species. The (p±) increases from pp and pe- 
ripheral Au+Au to central Au+Au collisions, especially 
for p, p and K . This behavior is consistent with an 
increase of radial flow with collision centrality. 

The K~ /ir~ and p/tt~ ratios of the integrated dN/dy 
yields within \y\ < 0.1 are depicted in Fig. 3(b). We 
observe little centrality dependence of the K~/tt~ or 
P/t:~ ratio from mid-central to central Au+Au collisions, 
indicating a similar freeze-out condition in these colli- 
sions. Similar centrality behavior has been observed for 
other particle ratios measured at y' snn = 200 GeV and 
130 GeV HE3. 

The observed centrality independence of K~ /tt~ is 
in contrast to low energy data at SPS 16] and AGS 
[l7j . where a continuous increase in K~ /it was observed, 
roughly doubling from peripheral to central collisions. To 
put our results into perspective with low energy data, 
we plot in Fig. 3(c) the K~/ir ratio as a function of 
(dN/dy) g ^ j n an a tt, em pt to reflect effects of both the col- 
lision energy and centrality. Here S is an estimate of 
the transverse overlap area based on the number of par- 
ticipants [|J, experimentally measured for the AGS and 
SPS data and calculated via the MC Glauber model for 
RHIC data [|| . The ( dN /dy)T ma y ^ e related to the initial 
conditions of the collision 0, E] , such as energy density. 



In high energy collisions the initial gluon density is sat- 
urated up to a momentum scale that is proportional to 

J ( rfjV Afyk 2] . Using data over a wide range of collision 
energy measured in various colliding systems, Fig. 3(c) 
shows a distinct change in the ratio behavior. Low en- 
ergy measurements (each representing approximately top 
60% of the geometrical cross section) appear to follow a 
trend that saturates at RHIC energies. One interpreta- 
tion of this is that strangeness production at low energies 
depends on how the collision was initially prepared, but 
not at RHIC energies. On the other hand, the K + /it and 
P/tt ratios do not reveal a common trend with ( dJV Afy)" _ 
However, we note that the net-baryon density, significant 
at low energies, greatly affects K + and p abundances 
through associated production of K + with baryons 
and baryon-antibaryon annihilation .19], respectively. 

In the framework of a chemical-equilibrium model 0, 
Efjj ]. integrated yield ratios can be described by a set of 
parameters: the chemical freeze-out temperature (T c h), 
the baryon and strangeness chemical potentials (/is, fi s ), 
and the strangeness suppression factor (7 S ). We fit our 
measured ratios with such a model to extract these pa- 
rameters. The value obtained for the chemical potential, 
fiB ~ 22 + 4 MeV, is independent of centrality within 
errors, and /i s is consistent with 0. The obtained 7 S in- 
creases from 0.56 ± 0.04 in pp to 0.86 ± 0.11 in central 
Au+Au collisions reflecting the measured K/tt ratios. 
The obtained T c h is summarized in Fig. 4 as a function of 
charged hadron multiplicity, together with T^m and (/?) 
extracted from the blast- wave model fit to our data. As 
seen in Fig. 4 ; y ^ djv A Mzi increases with centrality, T c h is 
independent of it, Tkin decreases and (/3) increases with 
centrality. This suggests that Au+Au collisions of differ- 
ent initial conditions always evolve to the same chemical 
freeze-out condition, and then cool down further to a ki- 
netic freeze-out dependent on centrality. The expansion 
of the system gives rise to collective flow. 

During expansion from chemical to kinetic freeze-out, 
entropy density drops approximately as T 3 |5l| , implying 
that the system size at kinetic freeze-out is at least a 
factor of of the size at chemical freeze-out. This 

l kin 

suggests a time span from chemical to kinetic freeze-out 
in central collisions is at least of the order of At ps — 
1)R/(3 S w 6 fm/c. Here we have taken R = 6 fm, the 
Au nuclei radius, as an estimate of the system size at 
chemical freeze-out. 

In summary, we have reported transverse mass and 
rapidity distributions of tt^ , K^, p and p for pp and 
Au+Au collisions at ^/snn = 200 GeV at RHIC. A 
boost-invariant region of at least Ay w 1 is developed at 
mid-rapidity for particles within our measured p± range. 
The spectra are well described by the blast-wave model, 
yielding a decreasing Tkin and increasing (/3) with cen- 
trality, reaching the values of Tkin =89+10 MeV and 
(/3) = 0.59 ± 0.05 in the 5% most central collisions. Par- 
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FIG. 3: (a) Mean transverse momentum of negative particles and (b) K /it and p/n ratios as function of the charged 
hadron multiplicity. Open symbols are for pp, and filled ones are for Au+Au data, (c) Mid-rapidity K~ /n ratio as function of 
(dN/dy)„ gy S t ema t j c errors are shown for STAR data, and statistical errors for other data. 
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(stars), T c h (circles) and Thin (trian- 
gles) and (b) (j3) as a function of the charged hadron multi- 
plicity. Errors are systematic. 



ticle ratios vary smoothly from pp to peripheral Au+Au 
and remain relatively constant from mid-central to cen- 
tral Au+Au collisions. The K~ /it ratio from various 
collisions over a wide range of energy reveals a distinct 
behavior in \ dN /dy)<r _ \ chemical equilibrium model fit 
to the ratios yields a T c h insensitive to centrality with 
a value of 157 ± 6 MeV for the 5% most central colli- 
sions. The drop in temperature from T c h to Tkin and 
the development of strong radial flow suggest a signifi- 
cant expansion and long duration from chemical to ki- 
netic freeze-out in central collisions. From these results 
the following picture seems to emerge at RHIC: collision 
systems with varying initial conditions always evolve to- 
wards the same chemical freeze-out condition followed 
by cooling and expansion of increasing magnitude with 
centrality. 
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